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ABSTRACT 


The diurnal variation of the thermosphere is the 
key feature for a theoretical understanding of the en- 
tire thermospheric structure and of its energy balance 
(heating processes, heat loss, conduction and convection). 
The earlier attempts (in 1962) in calculating the diurnal 
variation failed to yield the diurnal maximum of density 
and temperature at the observed time of 1400 hours. In 
these attempts no horizontal interaction was included. 

The recent calculations of the horizontal wind system, 
which is set up by the pressure gradients of the atmo- 
spheric bulge, showed that winds in the order of 100 to 
200 m/s in the altitude range from 200 to 400 Km can pro- 
vide a large horizontal energy transport. It cannot be 
decided yet whether the horizontal advection alone is 
sufficient to remove the whole discrepancy between the 
calculated time of the diurnal maximum (1700 hours local 
time) and the observed time (1400 hours) . 

Recent considerations on anomalous heat conduction 
through the magnetosphere have some bearing on the 
question whether the upper thermosphere is fully iso- 
thermal. This is important for the calculations of the 
diurnal variation since even a small temperature gradient 
at the thermopause can have a non-negligible effect, in 
particular, if the temperature gradient varies with 


local time. 



1 . Introduction 


In the course of the past ten years a detailed 
picture of the structure of the thermosphere and lower 
exosphere has emerged from satellite and rocket meas- 
urements. In particular, the vast number of measure- 
ments of the drag exerted by the atmospheric atoms and 
molecules on the satellite has revealed five different 
effects which cause variations of density and tempera- 
ture of the neutral gas. Furthermore, the data have 
provided some insight into the dependence of atmospheric 
properties on latitude and season. In this brief re- 
view we shall restrict ourselves to a height range of 
100 Km to 1200 Km, where most of the satellite-and 
high-altitude-rocket-measurements were made. 

The 5 effects are: 

1) the diurnal variation ("atmospheric bulge") 

2) the solar activity effect (27-day-variation) 

3) the solar cycle effect (11-year-variation) 

4) the geomagnetic activity effect 

5) the semiannual variation. 

For detailed information we refer to a recent review by 
W. Priester, M. Roemer and H. Volland (1967) and in 
particular, for new results on the seasonal-latitudinal 
variations to papers by L. G. Jacchia and J. Slowey (1967) 
and G. M. Keating and E. F. Prior (1967) . Furthermore 



a detailed review on the observational results of satel- 
lite drag analysis was given by L. G. Jacchia at the 
IQSY/COSPAR Symposium 1967 at London (to be published in 
the Proceedings of the IQSY. In this paper we shall 
only give a few illustrations to some of the aforemen- 
tioned effects. A satellite which nicely displayed the 
diurnal variation and the geomagnetic activity effect 
superposed on it was Injun III. L. G. Jacchia and 
J. Slowey (1963) determined from the drag-induced ac- 
celeration of this satellite the atmospheric density 
at an altitude of 250 Km (corresponding to the perigee 
altitude) for the time from December 15, 1962 through 
June 29, 1963. These data are given in the lower part 
of Fig. 1. During the observational period the local 
time of the satellite's perigee passed twice through 
the 24-hours. The local time is given as the lower 
x-axis of Fig. 1. On the top the dates of observation 
are given in the usual scale and in modified Julian 
dates (MJD) . The solid and the dotted curves represent 
the corresponding theoretical density values (from 
Harris and Priester 1962b) for the 1963 prevailing level 
of solar activity marked S = 90 and S = 100. The histo- 
grams in the upper part give the daily geomagnetic in- 
dices A and the solar 10.7 cm flux F. The black dots 
P 

represent the densities during geomagnetically quiet 
days (Ap 2) . These data should be compared with the 
theoretical curves. There are, however, no corrections 


made for either the correlation of the density with the 
10.7 cm flux ("solar activity effect") nor for the semi- 
annual variation. The latter shows up as low density 
values in December 1962, January 1963 and in June 1963. 
The most remarkable feature of Fig. 1, however, is the 
strong response of the density to geomagnetic storms. 

The quiet-day diurnal density variation of 250 Km 
in 1963 has an amplitude of a factor of 1.9. For years 
of high solar activity, for instance in 1958/1959 the 
amplitude was much smaller. This is borne out in Fig. 2, 
which gives the density distribution in the height range 
from 200 to 700 Km for diurnal maximum (1400 hours local 
time) and diurnal minimum (0400 hours local time) for 
two levels of solar activity (Fig. 2 is taken from 
CIRA 1965). The high activity is represented by a 10.7 
cm flux F = 200, the low activity by F = 75. It is re- 
markable, that the entire variation of the density at a 
height of 600 Km exceeds two orders of magnitude during 
the eleven-year solar cycle. 

As an example of the semi-annual variation we used 
the recent results obtained by G. E. Cook and D. W. Scott 
(1966) from the satellite Echo 2 during 1964 and 1965. 
These data show how strong the effect is at an altitude 
of 1130 Km for the rather low level of solar activity in 
1964, 1965 (Fig. 3). It also shows that the June, July 
minimum is lower, as usually observed, as the January 



minimum. In 1965 the summer minimum is misplaced due to a 
rather strong continuous increase of the solar activity 
level during June, July 1965 which masks the semi-annual 
minimum. There is no satisfactory theoretical explanation 
yet for the semi-annual variation. 

The recent investigations of high-inclination satel- 
lites have revealed more information on the dependence of 
density and temperature on latitude and season, which 
proved that the bulge center for altitudes below 600 Km 
follows the seasonal changes of the declination of the sun 
(Jacchia and Slowey (1967) ) . For altitudes above 600 Km, 
where helium is supposedly the dominant constituent of 
the neutral gas for periods of moderately low solar activ- 
ity, the density was found to be largest on the winter 
hemisphere at high latitudes. This was termed the 
"winter helium bulge" by Keating and Prior (1967) . The 
finding of this surprising behavior of the density was also 
confirmed by Jacchia and Slowey (1967) . 



2 . The structure of the thermosphere 

The basic feature of the thermosphere is the very 
steep increase of the temperature with height in the 
altitude range from 100 to 200 Km and the leveling- 
off of the temperature curve in the upper thermosphere 
(200 to 500 Km) . 

The steep temperature gradient is caused by the 

heating derived from the absorption of the solar XUV 

radiation, in particular, in the wavelength range from 
o 

100 to 912 A which was repeatedly measured by H. E. 
Hinteregger and his collaborators (see: Hinteregger et 
al.. Space Rec. V 1965). This part of the solar spec- 
trum originates primarily (or at least to 50 percent) 
in the solar corona-condensations above sunspot areas. 
This became evident when the early observations of 
satellite orbital periods (Jacchia 1958) revealed time- 
variations of the atmospheric drag, indicating variations 
in density and temperature in the thermosphere, and when 
it was found that these variations correlate with the 
solar decimeter flux (Priester 1959) . From the theory by 
M. Waldmeier and H. Mueller (1950) it was known that the 
so-called slowly varying component of the solar decimeter 
radiation originates in the coronal condensations. The 
correlation with thermospheric parameters then revealed 
immediately that the XUV radiation also mostly originates 


in the condensations, due to the rather close correlation 
to be expected on theoretical grounds between the deci- 
meter and the XUV flux. It furthermore proved that the 
XUV absorption is responsible for most of the heating of 
the thermosphere. Later in 1962 the observations with 
0S0 I confirmed that the integrated XUV flux and the 10.7 
cm radiation were proportional to each other. 

For altitudes above 200 Km the increasing mean free 
paths of the neutral air molecules and atoms cause a 
rapid increase of the heat conductivity which in turn 
must lead to an isothermal structure in the upper thermo- 
sphere provided that there is no appreciable heat influx 
from above, for instance by heat conduction through the 
magnetosphere. This general structure was pointed out as 
early as 1951 by D. R. Bates (1951) (see also M. Nicolet 
(1961), F. S. Johnson (1958)). 

That the temperature of the isothermal part of the 
upper thermosphere undergoes a large diurnal variation 
of several hundred degrees Kelvin came to light only when 
the large diurnal variation in density was discovered 
from the analysis of drag data by L. G. Jacchia (1959), 

S. P. Wyatt (1959) and W. Priester and H. A. Martin (1960) 
independently of each other . 



3) The diurnal variation of the thermosphere 


The diurnal variation of density and temperature has been 
observed in the height range from 200 to 800 Rn showing a 
steep increase in the morning, a maximum at 1400 hours local 
time, followed by a steep decline during the afternoon and 
a considerably less steep decline during the night with a 
minimum at about 0400 hours local time. These observations 
pertain to equatorial and moderate latitudes. 

The effect is caused by the heating due to absorption of 
solar energy and by the heat conduction of the neutral gas 
which conducts the energy down into the meso^pause (85 Rn) 
where it will be mostly eliminated from the atmospheric 
energy balance by reradiation processes. 

Attempts have been made to calculate the diurnal behavior 
of the thermosphere by Harris and Priester (1962a, 1962b, 
1963a, 1963b, 1965) and in a slightly different way by Ma- 
honey (1966) . 

A completely satisfactory solution can only be expected 
if one solves the following set of hydrodynamic and thermo- 
dynamic equations for the spherical earth: 



(conservation of mass) , 


?-jjr + y 9t * p 


- div^KgradTj + p cfiv v = <j[ 


Div 6T -t — 0 

(conservation of momentum) , 


( 1 ) 


at 


P = I ? T 

r M > 


(conservation of energy) , 


(ideal gas equation) , 


with the five unknowns density (0 , pressure p, temperature 
T and the components u, v, and w of the wind velocity 
v = v(u,v,w). g is the gravitational acceleration, © is 
the viscous stress tensor, c^ is the specific heat at constant 
volume, K is the coefficient of heat conductivity and q con- 
tains the solar heat flux and the heat loss due to infrared 


reradiation of the atmosphere. 
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Due to limitations of computer capacity and unknown boun- 
dary conditions, in particular of the wind systems, Harris and 
Priester (1962), following an earlier procedure by R. Jastrow 
and L. Kyle reduced the system in such a way that it depends 
on height and time only for the latitudinal range of the equa- 
tor. This implied the assumption that horizontal energy trans- 
port by winds was of minor importance. Nowadays this assump- 
tion can be no longer maintained. Nevertheless the simplified 
calculations provided useful insight into the energy balance 
of the thermosphere. 

The calculations were based on the following premises: 

1) Hydrostatic equilibrium is reached sufficiently rapidly 
at all heights and times. 

2) Diffusive equilibrium holds in the entire range from 
120 to 1000 Kn. 

3) The boundary conditions at 120 Mn (number densities 
of N , O , O, A, He, H and the temperature) do not vary with 
time. This condition, however, was abandoned in later versions 
of the investigation (Harris and Priester 1965, and as yet 
unpublished calculations in 1966 and 1967) in order to stu- 
dy the effects of time-dependent boundary conditions. 



4) Heat flow was permitted only in the vertical direction 
- no lateral conduction or advection. This condition, again, 
was eased, when the effects of lateral conduction and convec- 
tion in the equatorial plane were investigated. 

As a heat source the absorption of the solar XUV-radiation 
in the range from 100 to 1000 A was used with average values 
for the absorption cross sections of N 2 » 0 2 and O. It was 
proved that taking average values for the cross sections over 
the entire wavelength range did not alter the height distribu- 
tion of heat in any significant way. This proof was done by 
splitting the total range into five wavelength ranges with 
the corresponding cross sections. Also the investigation by 
Mahoney (1966) who used 32 ranges of wavelength and the cor- 
responding cross sections showed no significant differences 
in the resulting atmospheric parameters as compared to the 
calculations with one average cross section for each of the 
three main constituents. 

The heat conversion coefficient was taken to be 37 per- 
cent and independent of height and the chemical composition 
of the atmosphere. In our 1965 calculations we also added 
the heat provided by absorption of the radiation in the 
Schumann-Rung e range (1300-1700A) . This, however, has only 
an effect on the magnitude of the temperature gradient at the 
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lower boundary (120 Km) , but does not significantly influence 
the diurnal behavior of the thermosphere. 

Furthermore, the heat loss by reradiation from oxygen 

atoms was taken into account. As D. R. Bates (1951) had 

shown, only the radiation emitted by the transition from the 
3 3 

to the level of atomic oxygen is important. In 

general the amount of energy radiated away by this process 
is less than 10 percent of the absorbed XUV flux. Only close 
to the lower boundary altitude (120 Km) does the emitted radiation 
become comparable to the absorbed amount of the XUV radiation. 

The solution of the time-dependent heat conduction equa- 
tion in connection with the formulae for diffusive equilib- 
rium and using the heat sources explained above did not yield 
a satisfacotry agreement between the observed and the calcu- 
lated densities and temperatures. For the latter we should 
bear in mind that the temperatures could not be observed 
directly in the upper thermosphere. But it is believed that 
the temperatures inferred from the observed densities are 
accurate to within a limit of about 20 percent. 

A more detailed mathematical description of this result 
and the comparison between observed and calculated atmosphe- 
ric parameters can be found in SPACE RESEARCH III (Harris 



and Priester, 1963a). An illustration of the discrepancy is 
given in Fig. 4. The curve T represents the "observed” tem- 
peratures, the curve T £uv represents the calculated diurnal 
temperature variation in the upper thermosphere. (The curve 
T c will be discussed in a later chapter of this paper.) The 
temperatures correspond to a rather high level of solar acti- 
vity, characterized by a 10.7 cm solar flux of F = 200 flux 
-22 2 

units (10 W/m Hz ) . The total range of solar activity 
within the eleven-year cycle ranges from above F = 250 down 
to F = 60. 

From Fig. 4 we conclude that within the frame of the 
outlined theory an agreement between observations and theory 
cannot be obtained because the calculated maximum occurs at 
about 1700 hours instead of 1400 where it is observed. Fur- 
thermore, the requirements for the efficiency of conversion 
of absorbed solar XUV radiation into heat were prohibitively 

large. The total amount of flux required for conversion into 

2 

heat was 2.5 erg/cm sec. This was the amount which had been 
measured by H. E. Hinteregger (1961). Thus apparently the 
requirement is a 100 percent conversion efficiency, which is 
unacceptable. Later on the situation regarding the absolute 
flux values eased somewhat due to improved measurements which 
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pertain to lower levels of solar activity. While the measured 
flux values remained practically the same, the amount of heat 
required in the theory decreased proportionally to the decrea- 
sing solar 10.7 cm flux. Since, however, the observed beha- 
vior of the upper atmosphere during the entire decreasing phase 
of solar activity from 1958 through 1963 rather conclusively 
requires that the XUV flux should vary in proportion to the 
10.7 cm flux it is believed that the earlier measurements of the flux 
yielded values considerably too small. We shall have to wait 
for the next solar maximum for new measurements of the XUV 
flux in order to reach final conclusions on the problem. 

In Fig. 5 the diurnal distribution of the heat sources 
integrated over the entire height range is given. These 
integrals correspond to the amount of solar flux available 
for the heating of the thermosphere. The curve labeled No. 

4 with a total flux of 2. 5 erg/cm 2 sec would yield the curve 
T_„.. in Fig. 4 . with the maximum late by 3 hours. The curve 
No. 3 is made in such a way as to yield agreement with the 

observed data (curve T in Fig. 4). in 1962 it was concluded 

M 

from the then available measurements that the amount of XUV 

2 

flux available for heating was only in the order of 1 erg/cm sec 



(under the assumption of a conversion efficiency of 37 per- 
cent). This is represented by the dotted curve No. 1. Since 
this flux was highly insufficient we concluded in 1962 that 
there should be a "second heat source”. The second heat 
source required for agreement between observation and theory 
is given by curve No. 2 in Fig. 5. The curves 1 and 2 together 
yield curve 3. In 1962 it was believed that the discrepancy in the 
energy required could be provided by the solar wind. This idea 

was stimulated by the fact that during geomagnetic storms a 
considerable heating of the thermosphere takes place, the 
energy of which derives from the solar wind. 

In the light of the situation with regard to the direct 
flux measurements it is now the preferred theory that there 
might be as much flux available for heating as indicated by 
curve 4 for a solar activity level of F = 200. In order to 
shift the energy into curve 3 a horizontal transport mechanism 
is required. This mechanism is sought in a global horizontal 
advective flow pattern, which is set up by the pressure gra- 
dients between the maximum and the minimum of the diurnal 
atmospheric bulge. The curves 3 and 4 are made in such a way 
that the areas underneath the curves are the same. Thus 
they both provide the same amount of total energy. 


4 • The effect of a global wind pattern on the diurnal 

variation “ “ ~ 

A few years ago the wind speeds which might occur in the 
thermosphere were usually grossly underestimated. Then, from 
artificial cloud release measurements and from theoretical 
considerations, it became more and more obvious that winds 
with speeds exceeding 100 m/s might occur in the thermosphere 
and that the diurnal bulge ought to generate a global wind 
system with speeds in the same range. 

Wind patterns have been derived from the pressure gradi- 
ents of the diurnal bulge by R. S. Lindzen (1966), J. E. 
Geissler (1966, 1967), H. Kohl and J. W. King (1966), H. Kohl 
(1967) and by I. Harris (1966). 

I. Harris also investigated the effect of the horizontal 
winds on the shape of the diurnal temperature and density 
distribution. Since the available computer capacity was 
insufficient for a full treatment of the winds on a spherical 
earth, he restricted the computations to the equatorial 
plane. But these calculations are the only ones from the 
references mentioned above which include non-linear inertial 
terms. Due to the restriction to the equatorial plane the 
resulting corrections on the diurnal variation can only be 
regarded as indicative of the order of magnitude of the 
effect at the different hours of the day. Nevertheless, 
it yields some valuable insight into the more comprehensive 
problem. 

Within the computations of the winds there remained 
one problem which is not yet fully explored. That is the 


question of the effectiveness of the ion drag which pre- 
vents the neutrals from flowing freely. In particular, 
in the equatorial plane the ions will be rather effective 
in reducing the wind velocity of the neutrals in longi- 
tudinal directions , since the ions can be regarded as 
essentially fixed on the magnetic field lines of the 
earth, which is basically perpendicular to the equatorial 
plane. In his calculations I. Harris used for the number 
densities of the ions the electron density profiles as 
given by S. Chandra (1963) for daytime and nighttime 
conditions. He further applied S. Chapman’s (1956 ) expression for 
the collision frequency between the neutrals and the ions 
and assumed an average of half the momentum transferred 
in each collision. Under these assumptions the calculated 
wind velocities in the equatorial plane remained rather 
small and therefore had no significant effect on the di- 
urnal temperature distribution. 

In order to see under what circumstances the wind 
system would be effective enough to provide the appropri- 
ate shift in the diurnal curve, the effect of the ion 
drag was reduced artificially in the computations. But 
only when the ion drag was reduced to 10 percent of its 
original value as represented by the above stated as- 
sumptions the wind speeds reached values as high as 
200 m/sec at an altitude of 300 Km. The thus obtained 
velocity distribution in the equatorial plane is given 



in Fig. 6. The basic feature is an East to West wind 
in the morning and a West to East wind in the late 
evening . 

The effect of this wind system on the diurnal 

temperature curve is shown by the solid curve labeled 

T in Fig. 4. We can see that the horizontal advection 
c 

has reduced the amplitude of the diurnal temperature 
variation considerably. It also shifted the entire 
curve to earlier hours of local time. There remained, 
however, a rather curious maximum close to 1600 hours. 
Since it cannot be expected to have any better agreement 
between "observed" and calculated temperatures due to 
the simplifications in the calculations, in particular 
the restriction to the equatorial plane, it must be con- 
cluded that a global wind system can provide the required 
shift of the atmospheric bulge center to a time of ap- 
proximately 1400 hours local time, provided that the ion 
drag does not effectively prevent high winds from occur- 
ring. Of course, the ion drag effect will be maximum in 
the equatorial plane and will be much less preventive in 
the other directions. 



5 . Limitations of the thermospheric models 

A remarkable feature of the satellite drag measure- 
ments is the large reliability of the correlation of at- 
mospheric parameters with solar and geomagnetic activity 
and the stability of the atmospheric bulge, the maximum 
of which remains close to 1400 hours local time through- 
out the solar cycle (so far established for the seven- 
year period from solar maximum through solar minimum only) . 
We have to bear in mind, however, that satellite drag 
results provide density data, which are averaged over a 
considerable range in time and distance. Therefore in 
any comparisons of the models' predictions with instan- 
taneous measurements some caution must be exercized. This 
is particularly suggested by the larger scatter of direct 
measurements and by the recent discovery of two different 
wave phenomena. One type was found in the lower thermo- 
sphere by K. Mauersberger et al. (1967) , the other kind 
in the upper thermosphere by G. Newton et al. (1967). 

Another important aspect of the theoretical calcu- 
lations of the thermospheric behavior is the constancy 
of the boundary conditions at both the lower and the 
upper height limit, usually taken at 120 Km and at about 
1000 Km. The choice of the height of the upper boundary 
is rather unimportant. It can be selected anywhere be- 
tween 500 and 1000 Km without any significant influence 
on the thermospheric structure. 


Time-dependent variations of the boundary conditions 
can, however, be important. Diurnal variations of the 
density and temperature at a height of 120 Km are be- 
lieved to remain rather small with a temperature change 
confined in a range of less than 10°K. If so, this will 
alter the numerical values of atmospheric model parameters, 
but is not apt to influence the concept of the physical 
processes involved. 

At the upper boundary usually only a value of the 
temperature gradient is assumed; mostly taken to be zero 
because of the supposedly very effective insulation of the 
thermosphere from the solar wind by the magnetosphere. 
Recently, however, it has become more and more apparent, 
that in the rather unstable plasma of the magnetosphere 
the ionic heat conduction perpendicular to the magnetic 
field lines, might be greatly enhanced at the expense of 
the conductivity parallel to the field lines, with the 
effect that the "perpendicular" conduction could reach 
nearly the same order of magnitude as the "parallel" 
conduction (T. Tsuda (1967) , W. Priester, M. Roemer and 
H. Volland (1967), H. G. Mayr and H. Volland (1967)). 

If this is so, it might distort the isothermal structure 
of the neutral gas in the lower exosphere, yielding a 
temperature gradient in the order of 1 °K/km at heights 


of about 800 Km. 



In this context it is worthwhile to look at the 
entire temperature structure and also at the deviations 
from thermal equilibrium which are quite significant, 
even at heights as low as 200 Km. (J. V. Evans (1967)). 

The recent measurements by the IMP II satellite by 
G. P. Serbu and E. J. Maier (1966a, b) on electron and 
ion temperatures in the outer magnetosphere revealed a 
steep increase of both temperatures with height. Fig. 7 
gives the results of these measurements together with the 
profiles of the neutral gas temperatures in the thermo- 
sphere. The latter are given for three different levels 
of solar activity. The diurnal average temperatures only 
are plotted. 

The solid curve in the upper right part of Fig. 7 has 
been calculated for ionic heat conduction perpendicular to 
the magnetic field, in a fully stable plasma, assuming a 

4 

heat reservoir at a height of 5.10 Km with a temperature 
4 o 

of 7.10 K. This is thought to be provided by the con- 
tinuous flow of the solar wind, the temperature of which 

4 o 

is in the order of 7.10 K. A sufficiently good agreement 

is obtained if a temperature gradient of 2.65 °K/km is 

4 

assumed at the upper boundary at a height of 5.10 km. 

Thus the observed ion temperatures could be considered to 
be compatible with this kind of heat conduction. The heat 
flow associated with it, however, would be extremely small, 
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in the order of 5.10~ 15 erg/cm 2 sec. This, indeed, would 
provide a perfect insulation of the thermosphere. 

Since the plasma of the magnetosphere, however, is 
far from being in a stable state, we can expect a heat 
flow which is many orders of magnitude larger than the 
figure given above. As H. G. Mayr and H. Volland (1967) 
argue, it might even reach an order of magnitude which 
is comparable with the conductivity parallel to the mag- 
netic field lines. Under those circumstances the cal- 
culations for a stable plasma would no longer have any 
significance. But also the insulation of the thermo- 
sphere would be considerably less effective. This might 
then lead to a non- zero temperature gradient in the lower 
exosphere and upper thermosphere. Since temperatures, so 
far, cannot be measured directly in the lower exosphere 
(500 to 1000 Km) it is very important to derive tempera- 
tures from height profiles of the different atomic con- 
stituents, in particular, of oxygen and helium in the 
height range up to 1000 Km. A further interesting 
feature of the temperature structure is the rather steep 
increase of the ion temperatures at heights above 400 Km 
as it was revealed by the Thomson scatter technique 
(J. V. Evans 1967). These data could nave an important bearing 
on whether the structure of the neutral gas in the lower exosphere 


is isothermal. Evans explains the ion and electron 
temperature profile by heating due to fast photoelectrons 
which are being produced in the thermosphere in the 
absorption process of solar XUV-radiation. The basic 
argument for the photoelectron heating are the observed 
seasonal variations in the ion and electron temperatures. 



6 . 


Conclusions 


In this paper we have attempted to give an account 
of our present understanding of the thermal structure of 
the upper atmosphere in the range from 120 Km to about 
1000 Km and of its key-feature, the diurnal variation. 

The basic problem how to explain the density and temp- 
erature maximum at 1400 hours local time is still un- 
settled; The open questions are: 

1. Is the solar XUV flux sufficient to provide the 

ent ^ re energy for the heating of the thermosphere, even 
for the highest level of solar activity? 

2. If so, can then the horizontal advection set up by 
the pressure gradients around the atmospheric bulge pro- 
vide enough horizontal energy transport in order to ac- 
count for the 1400 hours maximum of the density in the 
upper thermosphere? 

3. How large is the heat flux which arrives at the 
upper level of the thermosphere (thermopause) from the 
solar wind through the magnetosphere under geomagnetically 
quiet conditions? Is the heat flux large enough to dis- 
tort the isothermal structure of the upper thermosphere 
and lower exosphere? If so, does the heat flux have a 
significant variation with the orientation of the earth r in 
particular with local time? 

These 3 complexes of questions must be answered 
before a sufficient understanding of the physical be- 
havior of the thermosphere can be reached. 
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FIGURE CAPTIONS 


Fig. 1 


Fig. 2 


Fig. 3 


z 


: Densities at a height of 250 km above sea level 

derived from Injun III by Jacchia and Slowey (1963) 
for the time interval from December 15, 1962 through 
June 1963. The data are plotted as a function of 
local time of the perigee. The corresponding dates 
are given at the top, also as MJD (Modified Julian 
Dates) . For comparison, the geomagnetic indices A 
and the solar 10.7 cm flux F are presented. The ^ 
numbers on the densities indicate the geographic 
latitudes of the perigree. The curves represent the 
Harris-Priester models S = 90 and S = 100. 


: Comparison of CIRA Model densities with data derived 

from satellite drag measurements (from CIRA, 1965) . 

The models are for high solar activity (Model 1 j_ 

F = 200) and very low solar activity (Model 2, F = 75). 
The observational data are the Bonn model 1961 for 
F = 200 and the data derived by Jacchia and Slowey 
for March 1963 corresponding to F= 78 and 04 hours 
local time from the 5 satellites Injun III, Explorer I, 
Explorer VIII, Vanguard II and Explorer IX. The 
maximum and minimum exospheric temperatures are stated 
on the right. 


: Variation of atmospheric density at a height of 1130 Km 

derived from Echo II by G. E. Cook and D. W. Scott 
(1966) . The densities display a pronounced semi-annual 
variation. 
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Fig . 4 


Fig. 5 


Fig . 6 


Fig . 7 


FIGURE CAPTIONS (CONT"D) 


: Diurnal variation of the temperature at the thermo- 

pause (500 km) . The abscissa is the local time t. 

The observed temperature T (dashed-dotted curve, 
inferred from density measurements) is compared 
with calculated temperature variations. The dashed 
curve (T euv ) is calculated under the assumption of 
heating Eythe absorption of solar XUV-radiation and 
no horizontal advection. The solid curve (T ) is 
obtained when the effect of horizontal advec€ion in 
the equatorial plane is included. 

: The heat source functions , integrated over the entire 

height range, as function of local time t. The 
integrals I correspond to that portion of the in- 
coming solar energy which is converted into heat in 
the thermosphere. The physical interpretation of 
the 4 curves is explained in the text. 

: The wind velocities in the equatorial plane for 4 

different heights from 150 to 400 km as function of 
local time t. These high velocities were obtained 
when the ion drag was considerably suppressed 
(see text) . 

: The temperature profile of the earth's outer atmo- 

sphere. In the height range z < 10 3 km the neutral 
gas temperatures (diurnal averages) are given for 
three levels of solar activity F = 275, 150 and 
65. In the upper right the electron temperatures 
(circles) and ion temperatures (circles with + signs 
inside) as measured by Serbu and Maier (1966a, b) 
are given. 
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